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Abstract 



^ . ABSTRACT In previous studies of the theory of Coulomb excitation, 

the term in the electromagnetic interaction which is quadratic in the 
vector potential has been ignored. In this paper, we use qualitative ar- 
guments and detailed numerical calculations to show that this quadratic 
term must be included when the bombarding energy is in excess of about 
Q ' 20 GeV per nucleon. 



1 Introduction 



^ \ The subject of relativistic Coulomb excitation has received extensive study in the two 

decades since its theoretical foundations were established in the classic work of Alder and 
Winther 0]. In recent years, this subject has been shown to be relevant to the practical 
question of the stability of beams of relativistic heavy ions, since some of the processes 
which lead to loss of beam ions are initiated by Coulomb excitation p, |^, ^. 

In the semi-quantal approach to relativistic nuclear Coulomb excitation, the relative 
motion of the projectile and target is treated classically. Indeed, it is usually assumed 
to be straight-line motion at constant speed v. The evolution of the internal degrees of 
freedom of each nucleus, under the influence of the classical electromagnetic field produced 
by the other nucleus, is then followed using quantum mechanics. The effect of this classical 
electromagnetic field on the nuclear charge and current densities has two components: one 
that is linear in the electromagnetic potentials, and one that is quadratic. All previous 
work on the theory of Coulomb excitation has been based on the linear term. However, it is 
known that there are some electromagnetic processes in which one cannot neglect the term 
quadratic in the potentials. Examples are Thomson scattering of photons by electrons [^, 
and the Zeeman effect in hydrogen atom states of high principal quantum number [H. 
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Our object in this paper is to determine whether the interaction term quadratic in the 
potentials must be included when we study the nuclear effects of the highly retarded 
electromagnetic field due to a rapidly passing projectile. 



2 The electromagnetic interaction 

We follow the usual approach of the semi-classical theory of Coulomb excitation, in which 
we solve the quantum-mechanical problem of the response of the target nucleus to the 
perturbation provided by the electromagnetic field of the passing projectile (ref.[|ll). In 
the absence of this perturbation, the target Hamiltonian has the form 

^0 = E ^p. ■ p. + wiO (1) 

3 

Here C, represents the internal degrees of freedom of the target. When the effect of the 
projectile electromagnetic field is included, the target Hamiltonian becomes 

H = E ^(P.-7A-(r;,t))2 + e,y.-Hr;,t) + l^(C) 

= + + (2) 



with 



^1 = -E^[p.-A^^Hr;,t) + A-(r;,t).p, 



+ e,¥.^^*(r;,t), (3) 
2mc2 



= E:^A-(r;,t).A-(r;,t). (4) 



Cj is the proton charge or depending upon whether the jth nucleon is a proton or neutron. 
Previous investigations (ref.0, |^) of relativistic Coulomb excitation have focussed only 
on the effect of the perturbation Vi, the term linear in Cj. This term is evaluated using 
either the Lienard-Wiechert potential for y9^*(r^-, t) and A^*(r^-, t), or the Fermi- Weissaker- 
Williams method of equivalent photons 10|. In both approaches, it is assumed that the 



perturbation produced by V2, is of secondary importance, and can be neglected. We note 
that the classical non-relativistic equations of motion, with the full Lorentz force, 

dv ^, , e ■, 
m— = eE(r, t) + -v x B(r, t) . 

dt c 

is obtained from the Hamiltonian (2.2) only if V2 is included. 

Our object in this study is to determine whether the V2 term in the interaction makes 
a numerically significant contribution to the Coulomb excitation amplitude. We will show 
in the following that for lead projectiles with kinetic energy below about 20 GeV per 
nucleon, the V2 term can safely be neglected. However, for kinetic energies at or above 80 
GeV per nucleon, the contributions from V2 are comparable to the those from Vi. 
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We assume that the perturbing electromagnetic field is produced by a spherically- 
symmetric projectile of charge Zpe moving along a trajectory given by 

r = b + t>tz, (5) 

where b is the impact parameter vector perpendicular to z, and v is the constant projectile 
speed. The scalar and vector potentials due to this projectile at the point r' of the target 
at time t can be taken to be the Lienard-Wiechert expressions (ref. 0) 



b - p'|2 + 72(t;t - z'Y 

A-\r',t) = %mr',t)z. (6) 
We focus our attention on the Fourier transform of the matrix elements of the perturbation, 

-e''^' <MVi{t) + V2{t)\4>^> 

-OO ih 

= + (7) 

Here (pa is an eigenstate of Hq corresponding to unperturbed eigenvalue The first- 
order Born approximation for the transition between (p^ and is expressed in terms this 
matrix element, with u) given its on-shell value of {Ej3 — Ea)/h. 



3 Orders of magnitude 

We begin with preliminary comparisons of the relative magnitudes of of matrix elements 
of Vi and V2 and of their dependences on bombarding energy and uj. 

Bombarding energy dependence. The Fourier transforms of Vi and V2 require, respec- 
tively, the following time integrals: 



Vi(ib-p'U') ^ I 



^ ^/lb-p'|2+72|i;t-2'|2 



V2(|b-p'|,^') = \%e^-'< 



I b — P' I 2 + 72 I _ 



/|2' 



The two t integrations yield 



Vi(|b-p'|,z') = iro(— |b-p'|), (8) 

nv 7^ 



V.(|b-p'|,/) = (9) 
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V^^ {uj) and V^"^ {u) are obtained by calculating the matrix elements of these time integrals 
between the nuclear states (pa and (pp. 

Ko(^\h — p'\) in Equation (3.1) diverges like log(7) as 7 cxd, whereas 7 e~~^^~^ ' 
diverges like 7. Thus as the bombarding energy becomes very large, the A ■ A matrix 
element grows relative to the matrix element linear in {(f,A). Alternatively, we can say 
that the extra factor of 

7 

» — p'p + 7^|ft — z'p 

in V2 increases the effect of retardation, which is most dramatic at large 7. 

UJ- dependence. For u — > 00, both expressions (3.1) and (3.2) decay exponentially. How- 
ever, since ii'o(x) ^J^e~^, there is an extra factor of 1/^^ in the fall-off of -K'o(:^|b—p'|) 

compared to e t"' ^ . Thus at large 7, where both matrix elements are appreciable, the 
A ■ A matrix element increases relative to the {(p, A) matrix element as uu increases. The 
extra retardation associated with the A ■ A term makes the electromagnetic pulse sharper, 
and so favors higher uj values in the Fourier transform of the matrix element. 

Relative magnitude. We can get a crude estimate of the relative importance of V2 and 
Vi by considering the ratio 



^A-(r,t)-A-(r,t)_,^(>-(r,t))^ e ^.t^^^ _ 



eip''''^{r,t) ev9''^*(r,t) 2mc'^ v w v ^ 

To make this a little more quantitative, let us give ip^^^{r, t) its maximum value of '-fZpc/b. 
Then the above ratio is 

1 

For Zp = 82 and b = lOFm, this becomes ~ .00637(1 — ^). This has a value of ~ .01 for 
Tp/A = 1 GeV (7 = 2.066), and ~ .34 for Tp/A = 50 GeV (7 = 54.3). Thus we expect 
that Vi can be neglected compared to Vi at projectile kinetic energies per nucleon of 1 
GeV or lower, but at 50 GeV, V2 may be of comparable importance to Vi. It will be seen 
below that this expectation is verified by detailed calculations. 

4 Matrix elements 

The evaluation of V^^ (uj) has been fully described in the hterature (ref . 0, ^ ) . We therefore 
turn our attention to Vp^{uj): 



V^'Jiuj) 



dte^-'i ^^S'i^'vt) ] ' PJ,M,-MM.W). (10) 



2%mc^ J- 
1 
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The transition charge density introduced in the last hne of Equation ( pUD can be expanded 
in terms of spherical harmonics of f': 



j 



leading to a multipole expansion of our matrix element: 



We disentangle the b and p' dependence of |b — p'| by using the expansion 

„tl^ (£-m)!!(£ + m)!! ^ ' 

which is valid when p' < b. To complete the evaluation of the r' integration in (0), 
we express pL(r')F^(f') in cylindrical coordinates by expanding p2,(r')Fj^(f') in terms of 
3-dimensional harmonic oscillator eigenfunctions 

PL(r')n^(f') = E r^Hr') (13) 

n 

and then performing a unitary transformation to products of 2-dimensional oscillator eigen- 
functions ipN^M^p', f') and one-dimensional oscillator eigenfunctions ipN^i^z'): 

The coefficients needed for this expansion are given in ref.([|ll|]). The p', 0', z' integrations 
can now be done in (0), with the final result 

,L-\M\_ 

^max 2 I 

xE(-l)"cH E {-IY^nX-)^pKM + p)\ <nLM\2p+\M\,2n+L-\M\-2p,M > 

n=0 p=0,l,2... ^ 

V V i-lVi^Y^'-' iM + i + 2j + 2n')\\ 

f^ok, ^S^v^^ j!(M + j)!(p-j)!(2ri')n(2£ + 2n' + l)!! • ^ ^ 
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In this equation we have set = 7r/2 (the projectile trajectory is in the y — z plane), and 
assumed that both u and M are positive. We get the remaining matrix elements using 
the symmetry relations 

u in is the harmonic oscillator size parameter h/{muJosc), which is used for all the 
oscillator eigenstates. It is chosen to give the best convergence for the expansion (|TBp of 
the transition charge density. 'iPnA^) is an oscillator eigenstate with argument ^, but with 
size parameter l/u. In the examples given below, we will consider Coulomb excitation of 
a giant quadrupole excitation in "^^Ca (ref.p). We use the Tassie model ref. |13|, which 
describes the 2"*" resonance as a one-quantum vibrational oscillation of an incompressible 
irrotational fluid. Explicit expressions for the transition charge and current densities are 
given in ref . . In this case, expansion of the transition charge density in terms of harmonic 
oscillator eigenstates is facilitated by the finite sum 



- nL+2K -!^^^L,'~,^ k\{2L + 2k + U/^ 

z/r')^ e ~Ytr(r') = — ' .1 ^ 



ft t — \\n 1 

('^ - '^)' 7n!(2L + 2n + l)!!2^+2-n 



5 Application to the exitation of vibrational states 

The V2 matrix elements calculated in Equation (4.5), added to the well-known matrix 
elements of Vi, can be used in a full coupled-channel calculation to yield a complete 
description of the Coulomb excitation process. Since our object here is only to assess the 
relative importance of V2, we will limit our attention to the exactly-solvable vibrational 
model. The main result of this model is that the probability, P{b,n), of exciting a state 
with n oscillator quanta when the impact parameter is b is given by a Poisson distribution 



P{b,n) 



where q{b) is the square of the on-shell matrix element of the interaction between the 
n = and n = 1 states. The Coulomb excitation cross-section for the population of the 
n— quantum state is then obtained from P{b, n) by an integration over b: 



/•oo 

cr„ = / 2Mb P{b,n) 

^ brrt i n 



The lower limit b^^^ is chosen to be large enough to ensure that only electromagnetic inter- 
actions contribute to the process. We apply this model to a hypothetical giant quadrupole 
oscillator band in "^^Ca, with huj=20 MeV. The B(E2, 2+) is assumed to have a 

value of 450 e^fm'*, which exhausts the energy-weighted sum rule. The projectile is ^^^Pb, 



and we have taken bmrn to be 12 fm 14 
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Figure 1 shows calculated excitation functions for populating the magnetic substates 
of the one-quantum 1 = 2 level, using Vi alone. The curves labelled by M = |1|, |2| 
correspond to the summed cross-sections for M = ±1, ± 2, respectively. The strong 
decrease of the M = cross-section at high bombarding energy is a result of cancellation 
of the contributions of the scalar and vector potentials to Vi of (2.3). At the highest 
energies, the M = ±1 cross-section is greatest. This is related to the applicability of the 
equivalent photon method in this region, since a shower of photons moving in the +z 
direction can only have M = ±1. 

Figure 2 shows calculated excitation functions, including the effects of both Vi and V2. 
Since the A ■ A interaction essentially involves the delivery of two photons to the target, 
arguments limited to one-photon exchange no longer apply. In particular, the M = 
cross-section is no longer suppressed at high bombarding energy. The large cross-section 
is due to the fact that the strongly-retarded potential can produce M = transitions at 
very large impact parameters. 

A polarization-insensitive measurement would yield the cross-section incoherently summed 
over M. This is plotted in Figure 3, for Vi alone, and for the full Vi + V2 interaction. It 
is seen that for bombarding energy per nucleon below about 20 GeV, V2 makes a rela- 
tively small contribution, but at 80 GeV the full Vi + V2 interaction yields about twice the 
cross-section of Vi alone. At higher bombarding energy, the V2 interaction dominates. 



6 Discussion 

We have seen that qualitative arguments and detailed numerical calculations strongly 
indicate the importance of the A ■ A contribution to the response of a target nucleus to 
the electromagnetic field of a highly-relativistic projectile. It should be emphasized that 
approximation methods such as the FWW method of virtual quanta provide convenient 
descriptions of the electromagnetic field, but do not, by themselves, include the effect of 
A ■ A which is quadratic in ej. 

An uncertainty of our analysis is the applicability of the non- relativistic Hamiltonian 
(2.2) when the external electromagnetic field is very strong, even when the protons in the 
target are moving slowly . If we attempt to derive this Hamiltonian by taking the non- 
relativistic limit of a proton Dirac equation, we need to assume that both the proton kinetic 



energy and e*/?^* are small compared to the proton rest energy (see, for example, [|T5|). But 
we have seen in Section III that the condition e^?^* << mc^ is violated for grazing collisions 
when the bombarding energy per nucleon exceeds about 50 GeV. Thus, if nucleons are 
to be regarded as slowly-moving Dirac particles, the discussion of their electromagnetic 
properties requires something more complicated than Equation (2.2) when we are dealing 
with the fields encountered in ultra-relativisitic collisions. 

The V2 contributions to the curves in Figures 2 and 3 suffer from a gauge ambiguity. 
The on-shell matrix elements of Vi are gauge invariant; the on-shell matrix elements of V2 
are not. Of course, if we use V2 in a full coupled-channel calculation, and arrive at exact 
Coulomb-excitation cross-sections, these will be manifestly gauge invariant. We note that 
a full coupled-channel calculation using Vi alone will not be gauge invariant, since both 
Vi and V2 are required in Equation (2.2) in order to yield a gauge-invariant Schrodinger 
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equation. The vibrational model used in this paper is not fully gauge invariant. The gauge 
we have used is that implied by the Lienard-Wiechert potential, and is widely employed 
in studies of Coulomb excitation. Although we cannot claim gauge invariance for our 
results, we believe that we have shown that, at bombarding energies in excess of about 
20 GeV per nucleon, it is important to include the effects of the A • A component of the 
electromagnetic interaction. 

References 

[1] K. Alder and A. Winther, Electromagnetic Excitation (North-Holland, Amsterdam, 
1975). 

[2] A.J. Baltz, M.J. Rhoades-Brown and J. Weneser, Phys. Rev. E54 4233 (1996). 

[3] S. Datz et al., Phys. Rev. Lett. 79 3355 (1997). 

[4] LA. Pshenichnov et al, CERN-preprint nucl-th/991208. 

[5] J.J. Sakurai, Advanced Quantum Mechanics (Addison- Wesley, Reading Ms., 1967) 

[6] S. Gasiorowicz, Quantum Physics, II ed. (John Wiley & Sons, New York 1996), Chap- 
ter 13 

[7] A. Winther and K. Alder, Nucl. Phys. A319 518 (1979). 
[8] B.F. Bayman and F. Zardi, Phys. Rev. C 59 2189 (1999) 
[9] J.D. Jackson, Classical Electrodynamics (Wiley, New York, 1975). 

[10] C. A. Bertulani and G. Baur, Phys. Rep. 163 299 (1988). 

[11] B.F. Bayman, M. Bernath and D. Bes, Phys. Rev. C41 581 (1990) 

[12] L. J. Tassie, Austral. Journ. Phys. 9 407 (1956). 

[13] T. Suzuki and D.J. Rowe, Nucl. Phys. A286, 307 (1977) 

[14] H. Emling, Part. Nucl. Phys. 33 792 (1994). 

[15] H. Bethe, Intermediate Quantum Mechanics, W.A. Benjamin, New York, 1964, p. 211. 
Figure captions 

Fig.l Coulomb excitation cross sections for one quadrupole phonon for the standard 
term Vi. The curves correspond to different values of magnetic quantum number transfer. 
More details can be found in the text. 
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Fig. 2 Coulomb excitation cross sections of one quadrupole phonon for the full inter- 
action Vi_ + V2. 

Fig. 3 Comparison of the cross sections for Vi and V2 summed over the magnetic 
quantum numbers. 
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